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Abstract 

We have carried out a numerical simulation of a domain- wall model in (2 + 1)- 
dimensions, in the presence of a dynamical gauge field only in an extra dimen- 
sion, corresponding to the weak coupling limit of a ( 2-dimensional ) physical 
gauge coupling. Using a quenched approximation we have investigated this 
model at l3s{= ^/g1) = 0.5 ( "symmetric" phase), 1.0, and 5.0 ("broken" 
phase), where gg is the gauge coupling constant of the extra dimension. We 
have found that there exists a critical value of a domain-wall mass rriQ which 
separates a region with a fermionic zero mode on the domain-wall from the 
one without it, in both symmetric and broken phases. This result suggests 
that the domain-wall method may work for the construction of lattice chiral 

gauge theories. 
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I. INTRODUCTION 



Construction of chiral gauge theories is one of the long-standing problems of lattice field 
theories. Due to the fermion doubling problems, a naively discretized lattice fermion field 
yields 2^^ fermion particles, half of one chirality and half of the other, so that the theory 
becomes non-chiral Several lattice approaches have been proposed, but so far none of 
them have been proven to work successfully. 

Kaplan has proposed a new construction of lattice chiral gauge theories via domain-wall 
models 0]. Starting from a vector-like gauge theory in 2k + 1-dimensions with a fermion 
mass term being a shape of a domain-wall in the (extra) 2k + 1-th dimension, he showed 
in the weak gauge coupling limit that a massless chiral state arises as a zero mode bound 
to the 2/c-dimensional domain-wall while all the doublers have large masses of the lattice 
cut-off scale. It has been also shown that the model works well for smooth back-ground 
gauge fields 0,^. 

Two simplified variants of the original Kaplan's domain-wall model have been proposed: 
an "overlap formula" and a "waveguide model" 0. Gauge fields appeared in these 
variants are 2/c-dimensional and are independent of the extra 2k + 1-th coordinate, while 
those in the original model are 2k + 1-dimensional and depend on the extra 2k + 1-th 
coordinate. These variants work successfully for smooth back-ground gauge fields as 
the original one does. Non-perturbative investigations for these variants seems easier than 
for the original model due to the simpler structure of gauge fields. 

However it has been reported that the waveguide model in the weak gauge coupling 
limit can not produce chiral zero modes needed to construct chiral gauge theories. In this 
limit, if gauge invariance were maintained, pure gauge field configurations equivalent to the 
unity by gauge transformation would dominate and gauge fields would become smooth. In 
the set-up of the waveguide model, however, 2/c-dimensional gauge fields are non-zero only 
in the layers near domain- wall( waveguide ), so that the gauge invariance is broken in the 
edge of the waveguide. Therefore, even in the weak gauge coupling limit, gauge fields are 
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no more smooth and becomes very "rough" , due to the gauge degrees of freedom appeared 
to be dynamical in this edge. As a result of the rough gauge dynamics, a new chiral zero 
mode with the opposite chirality to the original zero mode on the domain-wall appears in 
the edge, so that the fermionic spectrum inside the waveguide becomes vector-like. It has 
been claimed |^ that this "rough gauge" problem also exists in the overlap formula since the 
gauge invariance is broken by the boundary condition at the infinity of the extra dimension 
||9|JlO[|. Furthermore an equivalence between the wave-guide model and the overlap formula 



has been pointed out for the special case dlTl. Although the claimed equivalence has been 



challenged in ref. it is still crucial for the success of the overlap formula to solve the 
"rough gauge" problem and to show the existence of a chiral zero mode in the weak gauge 
coupling limit. 

How about original Kaplan's model ? In this model there are two inverse gauge coupling 
(3 = Ijg^ and [3s = ^/gl, where g is the coupling constant in (physical) 2 fc- dimensions and 
gs is the one in the (extra) 2k + 1-th dimension. Very little are known about this model 
except /?s = case where the spectrum seems vector-like. In the weak coupling limit, 

corresponding to the — limit in this model, all gauge fields in the physical dimensions 
can be gauged away, while the gauge field in the extra dimension is still dynamical and 
its dynamics is controlled by [3s. Instead of the gauge degrees of freedom in the edge of 
the wave-guide, 2k + 1-th component of gauge fields represent roughness of 2k dimensional 
gauge fields. An important question is whether the chiral zero mode on the domain-wall 
survives in the presence of this rough dynamics. The dynamics of the gauge field in this 
limit is equivalent to 2k dimensional scalar model with Ls independent copies where Ls is 
the number of sites in the extra dimension. In general at large [3s such a system is in a 
"broken phase" where some global symmetry is spontaneously broken, while at small [3s the 
system is in a "symmetric" phase. Therefore there exists a critical point and it is likely 
that the phase transition at [3s = Ps is continuous(second or higher order). The "gauge 
field" becomes rougher and rougher at smaller [3s- Indeed we know that the zero modes 
disappears at /?s = 0], while the zero mode exists at /5s = oo ( free case ). So far we 



do not know the fate of the chiral zero mode in the intermediate range of the couphng (3s- 
There are the following 3 possibihties: (a) The chiral zero mode always exists except Ps ~ 0. 
In this case we may likely construct a lattice chiral gauge theory in both broken [j3s > (3f) 
and symmetric(/3s < (5f), and the continuum limits may be taken at (5s — (^g- This is the 
beset case for the domain-wall model, (b) The chiral zero mode exists only in the broken 
phase (/3s > (3^). In this case we may construct a lattice chiral gauge theory only in the 
broken phase via the domain-wall method. This is unsatisfactory, since the chiral gauge 
theory in the symmetric phase, which is an important theoretical basis for various models, 
can not be described via the domain-wall method, (c) No chiral zero mode survives except 
(3s — oo. The original Kaplan's model can not describe lattice chiral gauge theories at all. 

It is very important to determine which possibility is indeed realized in the domain- 
wall model. Therefore, in this paper , in order to know the fate of the chiral zero mode, 
we have carried out a numerical simulation of a domain-wall model in (2 -|- l)-dimension 
with a quenched ^7(1) gauge field in the P — oo limit. Strictly speaking, there is no order 
parameter in a 2 dimensional U{1) model( XY model ). On a large but finite lattice, 
however, a behavior of the 2 dimensional model is similar to the one of a 4 dimensional 
scalar model. Thus, we hopefully think that useful informations about the fate of the zero 
mode can be obtained from such a toy model in (2-|-l)-dimensions. In Sec. 2 , we have 
defined our domain-wall model with dynamical gauge fields. We have calculated a fermion 
propagator by using a kind of mean-field approximation, to show that there is a critical 
value of the domain- wall mass parameter above which the zero mode exist. The value of 
the critical mass may depend on (3s, which controls the dynamics of the gauge field. In 
Sec. 3 , we have calculated the fermion spectrum numerically using quenched approximation 
at (3s = 0.5, 1.0.5.0 and at various values of domain-wall masses. We have found that at 
any value of three (3s there always exists the range of domain- wall mass parameter in which 
the chiral zero mode survives on the domain-wall. Our conclusion and some discussions are 
given in Sec. 4. 
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II. DOMAIN- WALL MODEL 



A. Definition of the model 



We consider a vector gauge theory in d = {2k + l)-dimension with a domain-wall mass 
term, which has a shape of a step function in the coordinate of an extra dimension. This 
domain- wall model is originally proposed by Kaplan p|, and a fermionic part of the action is 
reformulated by Narayanan- Neuberger ||^, in terms of a 2 /c- dimensional theory. The model 
is defined by the action 



S — Sn -\- St. 



(1) 



where Sq is the action of a dynamical gauge field , Sp is the fermionic action. Sq is given 
by 



n,ii>u s 

+ /?,5:5:ReTr[?7^,(n,s)], 

n,/! s 



(2) 



where /i, i> run from 1 to 2fc , n is a 2/c-dimensional lattice point , and s is a coordinate 
of an extra dimension. Ufj,i^{n, s) is a 2/c-dimensional plaquette and U^din^ s) is a plaquette 
containing two link variables in the extra direction. (3 is the inverse gauge coupling for the 
plaquette U^^ and (3s is the one for the plaquette U^d ■ In general , [3 ^ (3s ■ The fermion 
action Sp on the Euclidean lattice, in terms of the 2/c-dimensional notation, is given by 

Sf = ):Y.H i's{n)-i^ \Us,f,{n)i)s{n + /i) - Ul^{n - - /i) 

S 



+ E E Mn) \MoPr + mJpJ Mn 



n s,t 



E E ^^(^) \Us,M^s{n + /i) + Ul^{n - fi)ilJs{n - f^) - 2?/',(n) 



(3) 



nfj, s 



where s,t are an extra coordinates , Pr/l = ^(1 ± 72fc+i) , 

{Mo)s,t = Us4{n)5s+i,t - ci{s)5s,t 
{Ml)s,t = UlUn)5s-i,t - a{s)5s,t. 



(4) 
(5) 



Here Us,^{n), Us,d{n) {d = 2k + 1) are link variables connecting a site (n, s) to (n + /i, s) or 
{n, s + 1), respectively, Because of a periodic boundary condition in the extra dimension , 
s, t run from —Ls to — 1 , and a{s) is given by 



a{s) = 1 — mo sign 
mo 



[s + i) sign(L, -s-^] 



1 - mo (-| < s < L5 - I) 

^ ^ (6) 

1 + mo {-Ls - I < s < -i) , 

where mo is the height of the domain-wall mass. It is easy to check that the above fermionic 
action is identical to the one in {2k + l)-dimensions, proposed by Kaplan [0,0 . 

In weak coupling limit of both (3 and jSg , this model becomes free theory and can be easily 
analyzed. In free theory at < mo < 1, it has been shown that a desired chiral zero mode 
appears on a domain-wall( s = plane ) without unwanted doublers. Due to the periodic 
boundary condition in the extra dimension, however , a zero mode of the opposite chirality 
to the one on the domain-wall appears on the ant i- domain- wall , s = — 1 . Overlap 
between two zero modes decreases exponentially at large Lg. A free fermion propagator is 
easily calculated and an effective action of a (2 + 1) -dimensional model including the gauge 
anomaly and the Chern-Simons term can be obtained for smooth background gauge fields 

i- 

Domain-wall models, however, have not been investigated yet non-perturbatively. Main 
question is whether the chiral zero mode survives in the presence of rough gauge fields 
mentioned in the introduction. To answer this question we will analyze the fate of the chiral 
zero mode in the weak coupling limit for (3. In this limit, the gauge field action Sq is reduced 
to 

SG = f3sT.Il I^eTr [v{n, s)V^{n + /i, s)] , (7) 

s n 

where the link variable Us,d{n) in the extra direction is regarded as a site variable V{n, s){= 
Us,d{n))- This action is identical to the one of a (d — l)-dimensional spin model and s is 
regarded as an independent flavor. The action eq.(|^) is invariant under 



V{n,s)-^gis)V{n,s)g^is + l) , (gis) e G), (8) 

where G is the gauge group of the original model. Therefore the total symmetry of the 
model is G"^^", where 2Ls, the size of the extra dimension, is regarded as the number of 
independent flavors. We use this (reduced) model for our numerical investigation. 



B. Mean field approximation for fermion propagators 



When the dynamical gauge fields are added even on the extra dimension only, it is 
difficult to calculate the fermion propagator analytically. Instead of calculating the fermion 
propagator exactly , we use a mean-field approximation to see an effect of the dynamical 
gauge field qualitatively. The mean-field approximation we adopt is that the link variables 
are replaced as 



V{n, s) 



(9) 



where z is a (n, s)-independent constant. From eq.(^ the fermion action in a (c? — 1) 
dimensional momentum space becomes 



s,t,p \ fJ, ' / 



(10) 



V(p) 



2 -s,, (Mt(^)).,= (Mt(^))., + ^5.„ 



'111 



where V(p) = Y^'^Ji 2(cos]9^ ~ !)• Following Ref. it is easy to obtain a mean field fermion 
propagator on a finite lattice with the periodic boundary condition: 

r 1 ~^ 

G{p)s,t 



-I s,t 

I (-^ E 7mP/. + m] GL{p)s,t] Pl + U-iT. 7mPm + ] GR{p)s,t \ Pr 



, (12) 



with p/j, — sin(p^). For large Lg where we neglect terms of 0{e with c> 0, Gl and Gr 



are given by 



5e-°+l^-*l + {Al - B) e-°+(*+*) + {Ar - B) e-"+(2i»-5-t)^ (^^ t > 0) 

^^^-a+s+a-t ^ ^^g-a+(L,-s)-a_(L,+t)^ (s > 0, ^ < 0) 

^^^a-s-a+t ^ ^^g-a_(L,+.)-a+(L,-t)^ (s < 0, t > 0) 

Ce-"- + - C) e"- (^+*) + {Ar - C) e""- (2^-=+^+*) , (s, i < 0) 



(14) 



5e-"+l^-*l + (^i? - B) e-"+(^+«+2) + (A^ - B) e-«+(2i^-*-*-2)^ (^^ ^ > _i) 

^^g-a+(s+l)+a_(t+l) _^ ^^g-a+(La-s-l)-a_(Ls+t+l)^ > ^ < _-|^^ 

yi^e"-(*+-^)~"+(*+-^^ + ALe~'^-^^''^^^^^^'^+^^''~^~^\ {s < — l,t > —1) 

Ce-^-\^-A + (Aj? - C) e"-(^+*+2) + (^^ - C) e-"-(2i»+^+*+2), (s, ^ < -i) 

(15) 



where 



V(p) 

a± = z{l — ^ mo) = zh±, 



a± — arccosh 

Al-- 
B = 



p2 + ^2 ^ ^1 



1 



1 



(16) 
(17) 
(18) 
(19) 



2a+ sinh ' 2a_ sinh «_ 

Behaviors of Ar, B and C as p ^ are similar to the ones in free theory: They have no 
singularity for all z A behavior of Ax, is, however, different: As p — > Ax, behaves as 

1 



Aj 



[(1 - m,f\ + 0{p' 



-, (0<mo<l-^), 



Ami - [(-^^ - 1) -m, 
ArnQZ^p^ 



2i2 



-, (1 — z < mo < 1). 



(20) 
(21) 



A critical value of the domain-wall mass that separates a region with a zero mode and a 
region without zero modes is mg = 1 — z. Since A^ term dominates for 1 — 2; < mo < 1 in 
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the Gl (eq.(0) ) and Gr ( eq.([T5|)), a right-handed zero mode appears in the s = plane 
, and a left-handed zero mode in the s = Lg — 1 plane. For < mo < 1 — z the right- and 
left-handed fermions are massive in all s planes. Since the terms of Al^Ar^B and C are 
almost same value in this region of mo, a translational invariant term dominates in Gl and 
Gr, so that the spectrum becomes vector-like. 

If 2; — > 1 , the model becomes free theory. The propagator obtained in this section agrees 
with the one obtained in Ref . . In the opposite limit that 2; — , since there is no hopping 
term to the neighboring layers, this model becomes the one analyzed in Ref. [0 in the case 



of the strong coupling limit /3s = , and in Ref. , in the case that z is identified to the 
vacuum expectation value of the link variables. This consideration suggests that the region 
where the zero modes exist become smaller and smaller as 2; (1 — 2 < mo < 1) approaches 
zero. What corresponds to 2 ? Boundary conditions 2 satisfies are 2 = 1 at /^^ = 00 and 
2 = at /?s = 0. The most naive candidate [0 



IS 



z={V{n,s)). (22) 

But this is not invariant under the symmetry (j^). The other choice invariant under (|^) is 

2^ = (TrRe{1/(n, s)V\n + /i, s)}). (23) 

If eq. (^) is true, zero modes disappears in the symmetric phase, where (V"(n, s)) = 0, 
while, for the case of eq. (p3|), the zero modes always exist in both phases, since 
(TrRe{V(n, s)V\n + /i, s)}) is insensitive to which phase we are in. 



III. NUMERICAL STUDY OF (2-M)-DIMENSIONAL U(l) MODEL 

A. Method of numerical calculations 

In this section we numerically study the domain-wall model in (2 + l)-dimension with 
a f/(l) dynamical gauge field in the extra dimension. As seen from eq.(^ , the gauge 
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field action can be identified with a 2- dimensional U{1) spin model (with 2Ls copies). In 
(2 + l)-dimension, 7- matrices are Pauli-matrices , ai , a2 , cr^. 

Our numerical simulation has been carried out by the quenched approximation. Config- 
urations of U{1) dynamical gauge field are generated and fermion propagators are calculated 
on the configurations. The obtained fermion propagators are gauge non-invariant in general 
under the symmetry (|^). The fermion propagator G{p)s,t becomes "invariant" if and only 
if s = t. Thus, we take the s — s layer as propagating plane(~ "physical space"), and 
investigate the behavior of the fermion propagator in this layer. 

To study the fermion spectrum, we assume a form of eq.(|l2D for the fermion propagator, 
from which we extract Gl and G^. We then obtain corresponding fermion masses from 
G^^(p) and G~^^{p) by fitting them linearly in p^, since, from eq.(|l3D : 

Gl^ = p"^ + M ^ mj , (p^O), (24) 
G-i = p2 + MM"^ ^ mj , (p^O). (25) 

We take the following setup for 2-dimensional momenta. A periodic boundary condition 
is taken for the 1st direction and the momentum in this direction is fixed on pi = 0. An 
anti-periodic boundary condition is taken for the 2nd-direction and the momentum in this 
direction is variable such as p2 = {2n + 1)tt/L , n = —L/2, ...,L/2 — 1.) 

If mj = , we conclude that there is a zero mode, and if mj 7^ , there is not. 

B. Simulation parameters 

Our simulation is performed in the quenched approximation on x 2Ls lattices with 
L = 16,24,32 and Ls = 16. The coordinate s in the extra dimension runs —16 < s < 15. 
Gauge configurations are generated by the 5-hit Metropolis algorithm at f3s = 0.5, 1.0, 5.0. 
For the thermalization first 1000 sweeps are discarded. 

The fermion propagators are calculated by the conjugate gradient method on 50 con- 
figurations separated by at least 20 sweeps, except at Ps =5.0 on a 32^ x 32 lattice where 
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the number of configurations are 11. We take the domain-wall mass mo = 0.7, 0.8, 0.9, 
0.99 at = 0.5, mo = 0.3, 0.4, 0.5, 0.6, 0.9 at (3^ = 1.0, and mo = 0.1, 0.2, 0.3 at = 
5.0 . The boundary conditions in 1st- and 3rd(extra)-directions are periodic and the one 
in 2nd-direction is anti-periodic. Wilson parameter r has been set to r = 1. The fermion 
propagators have been investigated at s = 0, 8, 15. These s are the layers where we put 
sources. The layer at s = is the domain-wall , at s = 15 , the anti-domain-wall , and at 
s = 8 , neither. Errors are all estimated with the jack-knife method. 

C. Quenched phase structure 

As explained before the gauge field action of our model is identical to that of the U{1) 
spin system in 2-dimensions. Therefore, there is a Kosterlitz-Thouless phase transition and 
this system does not have an order parameter on the infinite lattice. On the finite volume, 
however, we take a vacuum expectation value of link variables as an order parameter using 
rotation technique: 

^=<|4E^(^.^)I>- (26) 

where L is the lattice size of the 1, 2-dimension. 

The defined vacuum expectation value v above is zero in the Kosterlitz-Thouless phase 
but f > in the spin-wave phase on the finite lattice. (Increasing the lattice size, however 
, decreasing the value of f . In the infinite lattice size, the value of v is zero for all gauge 
coupling.) Since we are interested in the dynamics of 4-dimensional theories, where the 
phase transition separates a symmetric phase from a broken phase , we have used this 2- 
dimensional system on large but finite volume as a toy model of 4-dimensional real world. 
Therefore, in this letter, we refer to the Kosterlitz-Thouless phase as the symmetric phase, 
and to the spin-wave phase as the broken phase. Fig. |I|(a) shows that, on a 16^ x 32 lattice, 
V behaves as if it was an order parameter. From Fig. |T](b) we consider that the system is in 
the symmetric phase at = 0.5, while in the broken phase at Ps = 1-0, 5.0 . 
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D. Fermion spectrum in the broken phase 



At (3s = 1.0 and 5.0, the system is in broken phase. Here we mainly discuss the result at 
/3s = 1.0 in detail. 

We first consider the fermion spectrum on the layer at s = 0. Fig. |^ is a plot of the term 
corresponding to — sin(p2) ■ Gl and — sin(p2) ■ Gji as a function of p2 at mo =0.3 and 0.5. ( 
Note we always set pi = 0.) This figure shows that, as p2 goes to zero, Gl seems to diverge 
at mo =0.5 but stay finite at mo = 0.3, while Gr stays finite at both mo. Next let us show 
Fig. ^, which is a plot of the Gj^^ and G]f^ as a function of p2 = sin^(p2) at mo =0.3 and 0.5. 
In the limit p2 — > 0, G]^^ remains non-zero at both mo, while G~[^ vanishes at mo = 0.5. We 
obtain the value of m'j, which can be regarded as the mass square in 2-dimensional world, 
by the linear fit in and plot mj as a function of mo in Fig. ^ The mass of right-handed 
fermion, obtained from Gj^^, becomes very small ( less than 0.1) at mo larger than 0.5, so 
we conclude that this critical value is mg ~ 0.5. Whenever the domain- wall mass is larger 
than this value , this model produces the right-handed chiral zero mode on the domain-wall 
at s = 0. 

On the anti-domain- wall {s = 15), on the other hand, the mass of left-handed fermion 
becomes less than 0.1 at mo larger than the critical mass mp ~ 0.5, as seen in Fig. ^ It is 
noted that chiralities between the zero modes on the domain-wall and the anti-domain-wall 
are opposite each other. 

Finally Fig. ^ shows that, on s = 8 , the layer in the middle between the domain- wall 
and the ant i- domain- wall , both right-handed and left-handed fermions stay heavy. 

A similar result at jSg = 5.0 on s = is given Fig. |^. 

From these results above, we conclude that the domain-wall model with the dynamical 
gauge field on the extra dimension {i.e. the weak coupling limit of the original Kaplan's 
model) can create the chiral zero mode on the domain-wall, at least in the broken phase. This 
suggests that the original Kaplan's model has a great chance to work for the construction 
of lattice chiral gauge theories in the broken phase. 
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E. Fermion spectrum in symmetric phase 



The system is in the symmetric phase at jSs = 0.5. The fermion propagator is analyzed in 
the same way as in the broken phase. However, for example on the s = layer, — sin(p2) - Gl 
and — sin(p2) ■ Gr show similar behaviors on a 16^ x 32 lattice, as seen in Fig. |[ Smaller 
lattice sizes, stronger the similarity, which makes analysis more difficult in the symmetric 
phase. To see a difference between the right-handed and left-handed fermions, we have to 
take larger lattice size such as L = 24, 32. 

In Fig. 1^, we have plotted mass mj of both modes at s = as a function of ttlq. Although 
a difference of masses between the right-handed and the left-handed fermions is very small, 
about 0.1 or less at itlq = 0.99, this difference stays finite as we increase the spatial lattice 
size L from 24 to 32. Therefore we conclude that the right-handed fermion becomes massless 
at niQ larger than 0.9, while the left-handed fermion stays massive at all tuq, so that the 
fermion spectrum on the domain-wall is chiral. 

In order to see that the difference of mass between the right and the left is really a signal, 
not a statistical fluctuation, we have plotted mj vs. mo in the case of putting a source at 
the anti-domain-wall s = 15 in Fig. |1^. We observe, at mo =0.99, a massless fermion of the 
opposite chirality to the s = zero mode and a finite difference of masses between the right 
and the left, which stays finite as we increase the spatial lattice size. 

Furthermore, in the case of s = 8, the right-handed fermion and the left-handed fermion 



stay massive at all mo, as seen in Fig. |11 



From these results above, as the same case in the broken phase, we conclude that the 
original Kaplan's model can create the chiral zero mode on the domain-wall even in the 
symmetric phase. 
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IV. CONCLUSIONS AND DISCUSSIONS 



Using the quenched approximation, we have performed the numerical study of the 
domain-wall model in (2+l)-dimensions with the U{1) dynamical gauge field on the ex- 
tra dimension. From this study we obtain the following results. There exists the critical 
value of the domain-wall mass separating the region with a chiral zero mode and the re- 
gion without it, both in the broken and the symmetric phases of the gauge field. At the 
domain-wall mass larger than its critical value a zero mode with one chirality exists on the 
domain-wall and a zero modes with opposite chirality on the anti-domain-wall, and none in 
the middle between the domain-wall and the ant i- domain- wall. 

These results strongly suggest that it is possible to construct lattice chiral gauge theories 
at all jSs except for jSg = via the domain- wall method, and continuum limits may possibly be 
taken at the critical value of jSg where the phase transition takes place. In (2 + l)-dimensions, 
however, the gauge field in /3 = oo limit is special since there is no order parameter and the 
phase transition is topological. Thus, to make a definite conclusion on the construction of 
lattice chiral gauge theories via the domain-wall method, we must study realistic (4 + 1)- 
dimensional model with U{1) or SU{N) gauge field in /3 = oo limit. Such models in /? = cxo 
limit have a phase transition characterized by an order parameter, a vacuum expectation 
value of the link variables in the extra dimension. 

Moreover, it is interesting and important to find an appropriate correspondence between 
the propagator obtained in the numerical simulation and the mean field propagator with 
tuned parameter z. So far it is not clear what physical quantity is corresponding to z. 
Since zero mode seems to exist even in the symmetric phase, the correspondence (|22D is 
qualitatively incorrect. On the other hand, we have found that the correspondence (^31) 
fails to reproduce rriQ quantitatively. Since mean-field approximations can not work well in 
the lower-dimensions in general, we must try to answer these questions studying (4 + 1)- 
dimensional U{1) or SU{N) models. 
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FIG. 1. (a) Vacuum expectation value of link variables v on a 16^ x 32 lattice as a function of 
Ps ■ (b) A volume dependence of the vacuum expectation values of link variables v as a function of 
1/L. 
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FIG. 2. — sin(p2) • [Gl]o.o a-i^d ~ sin(p2) • [Gi?]o,o in the fermion propagator as a function of p2 
with pi = at Ps = 1.0 on a 24^ x 32 lattice, for mo=0.5(open circles) and 0.3(solid diamonds). 



18 



4.0 



3.0 



2.0 



1.0 



0.0 



0.0 



I/GlCPs) vs. sin P2 

P3=1.0, 24'x32, 50 conf., s=0 



♦ 

© 



0.2 



0.4 



0.6 



sin^P2 



Fig. 3(a) 



♦ m„=0.3 
O m„=0.5 



o 



0.8 



1.0 



4.0 



3.0 



2.0 



1.0 



0.0 



0.0 



1/Gr(P2) vs. sin P2 

P3=1.0, 24'x32, 50 conf., s=0 



0.2 



0.4 0.6 
sin^p^ 

Fig. 3(b) 



O m|3=0.5 



0.8 



1.0 



FIG. 3. [Gl]o}) and [Gr\qJ^ as a function of sin^(p2) with = at /Sg = 1.0 on a 24? x 32 
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FIG. 4. ruf vs. mo at Ps = 1.0 on a 24^ x 32 lattice, in the case of putting a source on the 
domain-wall s = 0, for the right-handed fermion(open circles) and the left-handed fermion(solid 
circles) . 
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FIG. 5. nif vs. mo at /?s = 1.0 on a 24^ x 32 lattice, in the case of putting a source 
on the anti-domain-wall s = 15, for the right-handed fermion(open circles) and the left-handed 
fermion (solid circles). 
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FIG. 6. ruf vs. mo at /3s = 1.0 on a 24^ x 32 lattice, in the case of putting a source on s 
the right-handed fermion(open circles) and the left-handed fermion(solid circles). 
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FIG. 7. m/ vs. mo at Ps = 5.0 on a 24^ x 32 lattice, in the case of putting a source on the 
domain-wall s = 0, for the right-handed fermion(open circles) and the left-handed fermion(solid 
circles) . 
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FIG. 8. — sin(p2) • [Gl]o,o (open circles) and — sin(p2) • [Gij]o,o (solid circles) as a function of 
P2 with pi = at /3s = 0.5 on a 16^ x 32 lattice. 
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FIG. 9. m/ vs. mo at Ps = 0.5 on x 32 lattices with L =24(circles) and 32(squares) in 
the case of putting a source on the domain-wall s = 0. Open symbol stands for the right-handed 
fermion and Solid symbol for the left-handed fermion. 
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FIG. 10. m/ vs. mo at Ps = 0-5 on x 32 lattices with L =24(circles) and 32(squares) in the 
case of putting a source on the anti-domain-wall s = 15. Open symbol stands for the right-handed 
fermion and Solid symbol for the left-handed fermion. 
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FIG. 11. ruf vs. mo at Ps = 0.5 on x 32 lattices with L =24(circles) and 32(squares) in 
the case of putting a source on s = 8. Open symbol stands for the right-handed fermion and Solid 
symbol for the left-handed fermion. 
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